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ABSTRACT 


The purpose of . this study is the development and experimental simulator 
evaluation of superimposed predictor s>'mbology in computer-generated pictorial 
displays. The display under investigation, is a tunnel display for the 
four-dimensional commercial aircraft approach-to-landing under instrument 
flight rules. It is investigated whether more complex predictive information 
such as a three-dimensional perspective vehicle symbol, predicting the future 
vehicle position as well as future vehicle attitude angles, contributes to a 
better system response, and suitable predictor laws for the predictor motions, 
are formulated. Methods for utiliiing the predictor symbol in controlling 
tl»c forwai'd velocity of the aircraft in four-dimensional approaches, are 
investigated. 

Ilie simulator tests show, that the complex perspective vehicle symbol 
yields improved damping in the lateral response as compared to a flat two- 
dimensional predictor cross, but yields generally larger vertical deviations. 
Methods of using the predictor symbol in controlling the torward velocity of 
the vehicle are shown to be effective. Tlic tunnel display with superimposed 
perspective vehicle symbol yields very satisfactory results and pilot 
acceptance in the lateral control but is found to be unsatisfactory in the 
vertical control, as a result of too large vertical path-angle deviations. 

Ilie research is carried out in the framework of the Terminal Configured 
Vehicle (TCV) program at Ungley Ucscarch Center, Hampton. Va., as part of an 
ongoing research on integrated advanced display concepts for the future com- 
mercial aircraft flight-deck. Tlio development and quantitative simulator 



evaluation of the predictive infonnation, has been carried out at the Flight 
Control Laboratory of the Technion - Haifa, Israel. In parallel, all 
developments are implemented in the tunnel display software package at Langley 
Research Center, for operation in the TCV research cockpit, and in a later 
Stage, for operation in the TCV B-737 research aircraft. At Langley Research 
Center, a qualitative simulator evaluation has been carried out. 
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I . IRTRODUCTIKi 


Computer-generated pictorial displays facilitate the integration of control 
information in a format, analog to the through-the-windshield visual field. 

The tunnel display, in which the three-dimensional approach path is displayed 
perspectively as a winding and descending ’‘tunnel -in-the-sky” is useful in 
particular for following complicated curved trajectories. In previous work 
[1], C1980), the tunnel display was successfully inrplemented to the steep 
and curved helicopter approach-to- landing under instrument flight rules. It 
is shown in Ref. [1] that pictorial displays, without further augmentation, 
yield impaired system damping due to the lack of peripheral visual cues. It 
is also shown, that superimposed predictor s>Tnbology furnishes the system 
with the necessary damping cues. 

^fhis research deals with the implementation of the tunnel display to the 
fixed-wing commercial aircraft approach-to-landing and with the exploration 
and development of more complex predictive information. In a curved approach, 
the trajectory curv^ature constitutes the main forcing function to the pilot/ 
vehicle system and is responsible for a large part of the pilot activity. 

Tne predictive information serves the pilot in coping with this forcing 
function. Two predictor laws are considered: (1) a non-linear, basic 

predictor law, based on a circular future vehicle path which is fomulated 
such that it enables following steady cur\^es of the trajectory' with zero steady 
state error in the lateral deviation and (2) a linear, more sophistiCdted full- 
state predictor law which provides a more accurate prediction but yields a 
steady-state error in steady curves, A linear/non-1 inear predictor law is 


/ 



ft 
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formulated which combines both the advantage of a zero steady-state error 
given by the basic predictor law, with the advantage of a more accurate 
prediction given by the full-state predictor law. 

In addition to the future vehicle position, displayed by a flat two- 
dimensional predictor cross, the three dimensional perspective vehicle symbol 
also displays the future vehicle attitude angles. Since the intercept angle 
between the vehicle and tunnel trajectory is equivalent with the rate of 
deviation from the trajectory, the display of the future attitude angle is 

expected to contribute to the system damping. 

In addition to providing the necessary damping cues, the perspective 
vehicle symbol is utilized in controlling the forward velocity of the aircraft 
in four-dimensional approaches. This is accomplished by using the changes in 
predictor distance, resulting from changes in forward velocity as a control 
cue. The advantage of using the perspective vehicle symbol for controlling 
the fonvard velocity is that all control information, necessary for l-ateral. 
vertical and velocity control remains concentrated in the central .area of the 

display. 



II. DESCRIPTIQ?! OF THE DISPLAY 


The basic tunnel display configuration for the 4-D approach- to-landing, 
is shown in Fig. 1. The winding and descending three-dimensional approach 
path is presented to the pilot as a tunnel-in-the-sky, which is inertially 
fixed in space. In order to follow the desired approach path, the pilot 
must keep the vehicle inside the tunnel. The image shown in Fig, 1 is 
analog to the "through-the-windshield” visual field and shows the horizon 
(a) and the tunnel image with cornerlines (b). The tunnel cross-section 
is constant and square and of 300 ft width and remains at all time upright 
with respect to inertial space and thus parallel to the horizon. Analog 
to the natural visual field, a left bank is visualized by a clockwise 
rotation of the image about the monitor center and a nose-up pitch motion 
is visualized by a vertical downwards displacement of the image, perpendicula 
to the horizon. 

Basically, the square tunnel elements are drawn at 200 ft Intervals. 
However, to prevent clutter, all elements are omitted in this configuration, 
and only the interconnecting comerlines are shown. These lines are of 
200 ft length and the points at whch these lines are connected, appear as 
^^l^^ht spots. UTiile moving througli the tunnel these spots highly contribute 
to the impression of forward motion. 

Superimposed on the tunnel image is the perspective vehicle symbol (c). 
Tile center of gravity of the vehicle symbol indicates the predicted vehicle 
location, T seconds in advance, and the angular orientation of the symbol 




indicates the predicted attitude angles of the vehicle. This predictor 
syirbol is located at a distance D ahead of the vehicle, where D is 
predicted from the actual vehicle velocity, vehicle states and control in- 
puts, T seconds in the future. The wing-span of the vehicle synhol is 
identical to the tunnel width, i.e. 300 ft., which is about three times as | 

large as the wing-span of the actual aircraft. The bars (d) are positioned 
on the vertical axis of the vehicle symbol, and serve as a vertical reference. 

The distance between these bars is identical to the height of the square, 

I 

i.e. 300 ft. In contrast to the tunnel Inage, the vehicle symbol is not | 

\ 

a "wire-frame" structure. A wire-frame syribol, in which all lines are i 

i 

visible, is ambiguous and can equally well be interpreted as pointing towards : 

the observer, as well as away from the observer. An unambiguous vehicle ■, 

symbol is obtained by removing the "hidden lines" from the vehicle symbol j 

i 

fuselage. Mi efficient technique for removing the hidden lines is given j 

in the appendix. j 

The four comer "tick-marks" (c) indicate a cross-section of the tunnel, ^ 

which moves along, ahead of the vehicle, at the same distance D as the I 

predictor syniol. The solid square (f) is a cross-section of the tunnel 
as well but is ixjsitioned at distance Do “ T/Vo ahead, where Vo 1® the 
desired velocity. Since the solid square corresponds to the desired velocity 
and the comer tick-marks to the actual velocity, the velocity of the vehicle 
in A-D approadics is controlled by matching the tick-marks to the solid square. 

In addition, an increase in velocity and thus an Increase in predictor distance, 
manifests as an apparent thrinking in size of the predictor symbol and a 
decrease in velocity manifests as an apparent growing in size of the predictor 


I 
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synbol. In order to augment the difference between the tick-marks and 
the solid square, the tick-marks are drawn at three times the intensity 
of the solid square. Furthermore, the tick-marks arc blinked at a 3 Hz 
frequency, when the velocity error exceeds a 5 ft/sec threshold. 

A digital read-out of the vehicle velocity Cg) is displayed in 
the bottom center of the display. In Fig. 2 a-c this display configuration 
is shown for Do =• 900, Do “ 1500 and Do “ 2000 ft. 


Tunnel with Two-Dlmensional_Predictor_Cross. 

The tunnel display with a flat two-dimensional predictor cross is 
shown in Fig. 3 a-c for three nominal predictor distances. The height 
and width of the cross are identical to the tunnel square, i.e. 300 ft. 

The center of the cross is identical with the center of gravity of the 
perspective vehicle symbol. However, in contrast to the perspective si*mbol, 
the predictor cross remains at all times upright with respect to the display 
and thus does not display the future attitude angles of the vehicle. 


Roll-Stabi^zed_Tunn^ 

In Fig. 4 the roll-stabilized tunnel display is shown. In contrast 
to the roll-version of the display, the horizon and tunnel image remain 
at all times level on the display and the roll-motion is visualized by 
banking the vehicle symbol. A disadvantage of this configuration Is that 
the vehicle syn^Dol displays the predicted bank-angle rather than the actual 
bank-angle. The bank-angle information is correct in the steady state only 
since then the predicted and actual bank-angle arc Identical. 




In Fig. 5 a configuration of the rolL-stabilized tunnel display is 
shown, in which the tunnel eleinents are banked in curves. The trajectory 
bank— angle corresponds to the bank- angle which is required in a coordinated 
turn at the desired vehicle velocity. Thus, the trajectory bank-angle pro- 
vides a bank— angle corunand. Both in level flight, as well as in a steady 

1 

) 

I turn, the actual bank— angle will be identical to the coniinanded trajectory 

I bank-angle and the wings of the perspective vehicle syiribol, will be parallel 

i . . 

f 

: to the base of the tunnel square. 

\ In Fig. 6 the roll-version of the banked tunnel is shown. In a steady 

i 

1 turn botli the wings of the vehicle sytnbol as well as the square will be 

parallel to the base of the irronitor. In this situation, the inclined 
horizon provides the only actual bank-angle infornation. 

In Fig. 7 the roll-version of the banked tunnel is shown for a tunnel 

i width of 450 ft. The perspective vehicle synbol and tunnel cross-section 

are enlarged accordingly. 

i 

I 



III. 


EXPERIKENTAL EVALUATION 


A. OBJECTIVES OF THE EXPERIMEIITAL PROGPAM 

The objectives of the experimental program are: 

Cl) To compare the performance of the tunnel display with a 

3-D perspective vehicle symbol with the one with a flat 2-D 
predictor cross. 

C2) To evaluate the use of the perspective vehicle s>Tnbol with 

velocity tick-marks for controlling the velocity of the aircraft. 

(3) To evaluate a more complex full-state predictor law. 

C4) To Investigate the effect of the tunnel width. 

(5) To evaluate the effect of displaying the commanded bank-angle by 
banking the tunnel elements in curves. 

(6) To compare the roll- version with roll-stabilized version. 

B. EXPERimilTAL SYSTEli 

The experinental program was carried out at the Flight- Control Laboratory 
of the Teclinion, Haifa, Israel. During the summer of 1981 all developments 
were inplenented in the tunnel soft^^are package at Langley Research Center, 
Hampton, Virginia, for operation in the germinal Configured V^^icle (TCV) 
research cockpit. At Langley Research Canter preliminary qualitative 
simulator tests were carried out. 

A functional diagram of the experimental system at the Technion is 
shoi^’n in Fig. 8. The vehicle response was comnuted in tv;o parts: Linear 

computations, such as velocities and angular rates in the body coordinate 
system, were carried out by analog simulation, and non-linear computations , 




such as transfo mat ions from body to intertlal system as well as the computa- 
tion of the vehicle path were carried out digitally. Two EAI-580 hybrid- 
analog coTi^uters served for simulating the linear part of the lateral and 
vertical dynamics, disturbance input shaping, as well as simulating engine 
dynamics, servo-actuators and stability augmentation system. Analog 
random processes for the distrubance inputs, lere generated by a Hewlett- 
Packard H0I-3722-A noise generator. The analog signals were converted to 
digital by a RTP 7431/30 analog- to-digitai converter. All the digital 
conqjutations were performed by a Data General Corporation Eclipse computer 
with a 16-bit word length and 128-K extended memory. Vehicle path compu- 
tations were performed in floating point \>rith 32-bit precision. The integ- 
ration timing was controlled by a real-time clock, and the analog system 
was fully slaved to the digital system. About 80% of the Eclipse Central 
processing Unit CCPU) capability was devoted to the generation of the 
images of the various display configurations. Optimised, special pro- 
gramming techniques were developed, based on efficient assembly written 
subroutines, using 16-bit fixed point arithmatics, for obtaining a sufficient- 
ly fast update rate. The generated images, digitally coded in a sequence 
of move /draw commands, were translated into analog voltages‘by a Hewlett 
Pacl^rd HP- 1350A graphics translator, for drawing the stroke-written image 
on a Hewlett Packard HP-1310A catliode ray tube, with electrostatic deflection 
system measuring 19 in, diagonally. Vehicle notions, thus presented to the 
pilot, were utilized in generating the control commands which, in their 
turn, were imparted to the analog computers. 

Control manipulators consisted of a tv/o-axis spring-loaded control stick 
and an unloaded throttle control lever. The range of the cwo-axls control 



tv 
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stick was from +2 to -2 cm and the maximum spring torque was 0.3 for both 
axes* A forward stick displacement created a pitch motion and a lateral stick 
displacement a roll motion. Rudder pedals were not present and turn coordina- 
tion was carried out by the stability augra^ station system. The range of the 
throttle control lever was from +8 to -8 cm. 

Recorded time-histories of vehicle motions, control commands and 
disturbances were recorded on-line in extended memory of the CPU, and memorized 
after each run on a 100 .megabyte magnetic disk. The memorized time-histories 
were retrieved off-line for further reduction. Condensed results, such as run 
averages and scores were tabulated and printed on a line-printer; time-historic ^ 
and graphs, were plotted on the CRT screen and photographed for documentation- 
Ail condensed results were permanently stored on magnetic tape. 

c. T2';C7I?TICN or TiiE ExrcKj:>xrrrs 

The experiments were concerned with the approach-to-landing in the range 
from 30,000 to 1,000 ft from the touch-down point. A plan view of the desired 
trajactory is shown in Fig. 9 and the vertical descent profile along the 
trajectory is shown in Fig. 10. Tlie commanded velocity Fq was set at 
243.6 ft/sec over the complete approach path. A description of the two types 
of experiments chat were conducted is given hereafter. 

1, Trajectory following in the presence of random lateral and vertical gust 
disturbances; Tlie subjects were instructed to nininizc the lateral and vertical 
deviations from the trajectory with minimum control effort. Each run started 
from initial location 7 with an initial lateral deviation of 200 ft to the left 
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of the trajectory, a vertical deviation of 100 ft below the trajectory and zero 
intercept angle. Thus, control action was required from the subject, immediately 
after starting the simulation run, to bring the vehicle back on the trajectory. 
The lateral gust disturbance components v and p and the longitudinal 
gust disturbance components u and w were generated by passing band-limited 
white noise with a band-limit of 1.5 Hz, through first-order shaping filters. 

The RMS value of was 13.5 ft/sec, of p^ 0.086 rad/sec, of was 9.0 ft/sec 
and of was 12.2 ft/sec, and the break frequency of all shaping filters was 

0.2 rad/sec. 

Each run lasted 128 seconds, during which the means and auto-covariances of 
deviations, state variables and control commands, were computed. 


2. Entering the trajectory from a random chosen, unknown location outside the 
trajectory: The experiment attempted to simulate a sudden confrontation with 

the situation of being located outside the trajectory. The subjects were 
instructed to bring the vehicle back on the desired trajectory, as fast as 
possible, as smooth as possible and with minimum control effort. Gust 
disturbances were not present in this experiment. In order to prevent the 
subject from knowing his initial position before the start of a simulation run, 
the display was initially blanked and was made visible only immediately after 
starting the simulation run. 

Each run started randomly from one of the 6 initial locations shown in 
Fig. 9. For all locations the initial lateral deviation was 2,000 ft to the 
left or to the righ of the trajectory and the initial vertical deviation was 
300 ft above or below the trajectory. The initial intercept angle was set 
between 0 and 60 degrees. 




\ 
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Each entry run lasted 38,4 seconds during which the following performance 


scores were computed* 1. The lateral settling time defined as the time 

from the start of the run to the moment the lateral deviation settles within 

a ±100 ft settling tolerance about the desired trajectory and the vertical 

settling time T , for which the settling tolerance is ±75 ft* 2* The means 
s V 

and auto-convariances of deviations, state variables and control commands. 

3. The lateral drviation score, defined as the averaged absolute value of the 
lateral deviation, computed between t«!Ti and t-Tzf where Ti « 10 sec and 
Tz «38*4 sec, according to: 

T2 

ti^Ti 

The vertical deviation score was computed in the same way as the lateral one. 
The lateral and vertical deviation scores were chosen to be averaged absolute 
values rather than averaged squared values, in order to prevent these scores 
from being dominated by the large initial deviation. For the same reason the 
averaging process only started at sec. Ti is chosen to be about 25% 

less than the best possible settling time. 

D. RESULTS 

Four subjects participated in the experimental program. Apart from 
subject A, all subjects were male. Subjects A, C and D were Aeronautical 
Engineering students with no prior flight or simulator training, and subject B 
an Aeronautical Engineer v;ith extensive simulator experience. Each subject 



participated in 2-3 2-hours /week simulator sessions. All subjects reached a 
stable well-trained level after 6 weeks of training. Subject motivation was 
largely enhanced by using a reward system, based on a general performance 
score, which was displayed to the subject at the end of each run. This score 
was composed of the weighted sura of mean-squared deviations and control 
commands. The reward, given in the form of extra pay, was determined by the 
performance level which was reached and maintained during the session. 

The results of the four subjects are treated separately and are 
su mma rized in Tables I-V. The results for each experimental condition 
represents the average and standard deviation of a set of six or more repeti- 
tions. 

D.l Results of Trajectory Following 

Comparison of the 3-D perspective vehicle symbol with flat 2-D 
predictor cross. 

The lateral results of the four subjects are shown in Figs. lJ-14. 
Both for the perspective vehicle symbol, as well as for the predictor cross, 
the lateral deviation was found to increase strongly with the predictor 
distance Z?, whereas the roll activity and lateral stick activity were found to 
decrease strongly. Tlie contribution of the perspective vehicle symbol was 
i‘ou;:d in particular in the significantly lower roll-activity, for all 
four subjects and almost over the complete range of D, which proves that the 
perspective vehicle symbol contributes to the system damping- 


On the other 


hand I the perspective vehicle cyinbol yielded a somewhat larger lateral 
deviation. In the vertical control theso cffecto wore less pronounced. T!ie 
vertical results of the four subjcctn arc shown i\\ Figa. 15-18 and indicate 
that also the vertical deviation strongly increased with IK The perspective 
vehicle s>Tnbol yielded a Koroevhat larger vertical deviation than the predictor 
cross, which might be attributed to the fact that the vehicle Bviabol with its 
complex shape is harder to match to the tunnel square. 

Time histories of single runs of subject B, for the perspective vehicles 
symbol and for the 2-D predictor cross, arc compared in Figs. 19-24, Figs. 

19 and 20 clearly demonstrate the lower roll activity as well aa lateral at ick 
activity for the perspective vehicle H>Tnbol at I’o- 1,500 and Po - 900 ft. At 
2,000 the difference is loss pronounced, probably because the lateral 
motion ot the vehicle s\Tnbol provides sulficionc damping cues, see Fig, 21, 

T!\e time-histories v the vertical control, shown in Figs. 22-24, 
indicate larger vertical deviations for the perspective vehicle s>Tibol, 
where. IS the pitch and stick, activity .are the same as for the predictor cross. 

- « L‘i i l2£iL'* c o n t r o 1 

The results of auto-throttle control and manual velocity control by means 
ot velocity tick-marks, are showTi in Figs, 11-18, Neither in the lateral, nor 
in the vertical control did tlie results for tlie runual throttle differ 
signif icant ly from the results for the auto-throttle, which proves that 
velocity control conld he carried out witliout affecting the tunnel fw'illowing 
perlox-:tiance. For comparison, manual velocity control was carried cut by means 





Results of the banked tunnel. 

Tlie banked tunnel was investigated, both in the roll-version, as veil as 
in the roll-stabilized version. Tlie performance of the roll-version of the 
banked tunnel was very similar to results of the straight tunnel, see Tables 
1-lV. Tlic subjects cornmented, that in the roll-version, the bank-angle conmiand 
information, provided by the tunnel elements banked in curves, did not really 
contribute to the following performance, and was ignored in most cases. 
Furthermore, the bank-angle command was only correct after entering a steady, 
coordinated turn. The incorrect bank-angle conimand in transients to cur\^cd 
sections was found confusing. 

The roll-stabilized version of the banked tunnel yielded generally larger 
lateral deviations and roll-activity, than the roll-version, see Tables I-IV. 
This was attributed to the fact that the actual bank-angle, which in the roll- 
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version is displayed by the inclination of the horizon, is no longer available 
in the roll-stabilized version. The predicted bank-angle, displayed by the 
vehicle syxjbol was clearly not sufficient. These results are confirmed by 
comparing the time-histories in Fig. 28 for the roll-stabilized banked tunnel, 
with the tine-histories in Fig. 20 and 23 for the roll-version of the straight 
tunnel. Subject opinion of the roll-stabilized version of the banked tunnel was 
less favourable than the one of the roll-version of the banked tunnel. 

The cf fe et of the . tunnel width . 

The results of the tunnel of 450 ft width varied between the subjects, see 
Fables I -IV. For subject B the roll-activity was significantly lower and the 
lateral deviation significantly larger for the 450 ft tunnel than for the 
300 ft tunnel. Also for subject A the roll-activity was lower for the 450 ft 
tunnel, but the lateral deviation was the same. However, both subjects A and 
B shoved a lover predicted lateral deviation for the 450 ft tunnel, which 
indicates that the subjects devoted more attention to the lateral error 
between vehicle s>Tnbol and tunnel square. Also subject C shoved a lower roll- 
activity for the 450 ft tunnel, and smaller lateral deviation. The subjects 
commented that generally the 450 ft tunnel was easier to control than the 
300 fc tunnel. 

full-state predictor law. 

The full-state predictor law yielded slightly larger roll-activities and 
sligh.tly Ivirger lateral deviations than the basic circular path law. Tliis 
might be attributed to the fact that in the full-state law the bank-angle rate 
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component of the bank-angle prediction, is much smaller than in the basic law, 
sec Eqs. (A. 96) and (A. 97) of the Appendix. However, subject opinion was 
more in favour of the full-state predictor law, since disturbingly rapid roll- 
motions found with the basic law, were missing. 


D.2 Results of Trajectory Entry 

Tlio entry experiment was conducted in a series of six runs. For each run 
the initial location was chosen at random and without replacement from the 
set of six initial locations given in Fig. 9. Significant differences between 
the scores were observed, for t!^c various initial locations. In order to rate 
the general entry performance, a series score was computed by averaging the 
results of the 6 runs in each series. For eacl^ display configuration at least 
6 series of runs were performed. The experimental results are surma ri zed in 
Table V and Figs. 29-32 and represent the average and standard deviation of 
sets of at least 6 series-scores . 

The roll-activity and lateral stick activity tended to be lower tor the 
nerspectivc vehicle s\m\bol than for the 2-D predictor cross. Note the 
Jiffcrence in control strategy between the subjects: in contrast to subjects 

V <ind n, subjects A and B show less control activity on the account of a 
larger lateral score, see Fig. 29. 

Tb.e vortical results in Fig. 30, clearly show a larger vertical score for 
the perspective vehicle symbol, which might be explained by the fact that the 
perspective s\Tibol is harder to match to the tunnel square. Time-histories of 
single entries arc shown in Figs. 31 and 32. No significant difference in conttv>l 
strategy between the perspective vehicle s\*mbol and the 2-1) predictor cross are 


not i cod , 




D.3 Results of Preliminary Simulator Tests at Langley Research Center 

In the qualitative evaluation at Langley Research Center the tunnel with 
perspective vehicle symbol and basic circular path predictor law, was investigated. 
In most runs, no atmospheric disturbances were present. The pilot subjects found 
the perspective vehicle symbol more difficult to familiarize with, than the flat 
2-D predictor cross. After familiarization with the display, the predicted pitch 
and yaw attitude angles, were found useful. However, the predicted bank-angle was 
found to be confusing since it was composed of both actual bank-angle as well as 
bank-angle rate. It was found preferable to bring the bank -angle rate portion of 
the prediction to zero and to display the actual bank -angle only. Note that for 
the full-state predictor law, the bank-angle rate portion of the bank-angle 
prediction is considerably smaller than for the basic predictor law, see Eqs . 

(A. 96) and (A. 97) of the Appendix. The pilot rejection of the bank-angle rate 
portion of the prediction, thus confirms that the first order bank -angle predic- 
tion of Eq. (A. 34) is not sufficient. 

The following parameters were varied: (1) Nominal predictor distance Do^ 

(2) Tunnel width, (3) Predictor size and (4) Trajectory curvature. At a nominal 
airspeed of 130 knots, the most suitable predictor distance was between 900 and 
1,250 ft. A tunnel width of 450 ft and a perspective vehicle s\Tnbol with a wing 
span of 80% of the tunnel width, was found adequate. 

Two trajectories were tested: a strongly curved path yielding steady state 

bank-angles of about 27.2 degrees and a moderately curved path, yielding steady 
state bank-angles of about 16.6 degrees. Both for the strongly curved path, as 
well as for the moderately curved path, the perspective vehicle symbol yielded an 
adequate perfornumce. However, the maximum bank-angle for the strongly curved 
path was too large to be acceptabJe in actual flight. 
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Manual velocity control by means of velocity tick-marks proved successful. 
Accurate velocity control was obtained without overshoots and with minimal 
throttle activity. The manual velocity control did neither significantly 
increase pilot workload, nor affect the tunnel following performance. 

Tlie banked tunnel was investigated, both for the roll-version, as well 
as for the roll-stabilized version. The roll-stabilized version was found 
preferable due to the fact that the bank— angle command was directly 
perceived as the inclination of the tunnel elements with respect to the 
stabilized horizon, and thus only with respect to the monitor frame as well. In 
contrast with the simulations at the Technion in which the predicted bank- 
angle was displayed, at the Langley version the actual bank angle was 
available, since the bank-angle rate portion of the prediction was set to 
r.ero. Favourable opinion was given to the fact that in a steady coordinated 
turn the wings of the perspective vehicle symbol were parallel to the tunnel 
square. 

TWO main problems were encountered: (1) The vertical' control was too 

inaccurate and yielded unacceptably large vertical path-angle variations; 

(2) In the lateral control transients to and from curved sections of the 
path uece coo sudden. 

The reasons for the unsatisfactory vertical control are two-fold: 

(a) Since a contact analog display is used, horizontal and vertical picture 
scales are identical, i.e., *45 degrees of visual angle. The visual angle 
which is satisfactory in the lateral control, is Coo insensitive in the 
vertical control. However, a reduced vertical visual angle is not possible 
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since it would seriously distort the image, (b) The path-angle is not 
explicitly displayed. 

Pilots commented 9 that although the lateral control was very satisfactory 
in straight sections as well as in steady curves, the transients from straight 
to curved sections were too sudden. This resulted from the fact that the 
curvature along the trajectory was varied in steps, without transients from 
one section to the other. 
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IV. COnCLUSlOHS AND SUGGESTIONS FOR FURTHER RESEARCH 

1. The perspective vehicle symbol requires a longer time to familiarize with 
than the 2-D predictor cross. 

2. The perspective vehicle symbol yields a smaller roll-activity and thus 
contributes to the lateral system damping. 

3. The perspective vehicle symbol yields generally larger vertical deviations. 

4. The perspective vehicle symbol with velocity tick-roarks enables accurate 
velocity control, without overshoots, with moderate throttle activity and 
without increasing the pilot workload or affecting the tunnel following 
performance • 

5. A tunnel width of 450 ft and perspective vehicle svTtbol with a wingspan of 
80% of the tunnel width yields the best pilot acceptance. 

6. The full-state predictor law yields a better pilot acceptance than the 
basic circular path law, in particular with regard to the bank-angle 
prediction. 

7. Tlie vertical control in the present display configuration is too inaccurate 
and yields too large vertical path-angle variations. Future research 
efforts should be directed to the augmentation of the display with super- 
imposed symbology, explicitly displaying the vertical path-angle in the 
correct scaling. The s^Tnbology should be located such, that it does 
neither interfere with the lateral control, nor clutter with existing 
symbology and yet is located as centrally as possible on the display. 

8. The banked tunnel in the present configuration, is only effective in a 
steady coordinated turn. Incorrect bank-angle corraands in transients are 
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found confusing. Future developments should deal with the definition 
3 higher— order continuous function for the trajectory curvature and 
commanded bank-angle. This function should be custom-taylored to the 
average vehicle response, in order to enable the pilot to follow the 
trajectory through straight, transient and curved sections by matching 
the actual bank-angle to the commanded bank-angle. 
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A. VEHICLE MODEL 


V. APPO^DIX 






The lateral and longitudinal dynamics are assumed to be fully decoupled. 
The linearized lateral equations of motion of the vehicle, for small devia- 
tions from the nominal trim condition, are given by: 



(A.l) 


where: y* = = y^/c/, 


and the “primed** derivatives Ll and Ni are given by 
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where i can be P. i\ 6 , 6 • 

^ a r 

In this study, the stability derivatives of a DC-S aircraft are chosen, 
in the approach-to-landing, trimmed at a nominal airspeed of 243.5 ft /sec and 
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with fl^ps pxtopded to 35 destoo,. The dimensional stability derivatives 
are derived from Ref. [ 2 ] and given in Table VI. 

The lateral stability is largely improved by a lateral stability augments 
tion system including a roll-rate on aeleron feedback, a yaw-rate on rudder 
feedback with wash-out and a sideslip on rudder feedback. Aeleron and rudder 
servo actuators are modelled by a first-order lag with a 10. rad/sec break 

frequency. Tlie equations of motion of the control system, including sen^o 
actuators, arc given by: 


6^ » 

“10.05 

+ 10.0 


a 

i? 

s 



0 =* 

>> 

“10.05 

>* 

+ 10.0 
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* s 



where and c. are the inputs to the sei-vo actuators given by;. 

“ l.O.v + 

^ G 

c = -8. Op + 8.0 >• 

w.ierc i;, the state of the wash-out filter, given by 

= -1.0i> + I* 


(A. 3) 
(A.4) 

(A. 5) 
(A. 6) 


(A. 7) 


and 0 ^, is the conmand input to the system, originating from a lateral control 
stick motion. Turn coordination is accomplished by means of the sideslip on 
rudder feedback and the rudder pedals are not used. 

The lincariced longitudinal equations of notion of the vehicle, for small 
deviations iron the nominal trim condition are given by: 




The dimensional stability derivatives are given in Table VI. 


I 

The longitudinal stability is augmented by a pitch rate on elevator feed- 1 

back and an angle-of-attack on elevator feedback. Tiie elevator servo is 
modelled as a first order lag and the equations of motion of the control system 
are given by: 

= -10.0 6^ + 10.0 e (A. 10) j 

where: = 6.0 q + 0.05 u + 6^ (A. 11) 

Throttle and engine dynamics are modelled by a first order lag with a 
time constant of 10 seconds. The auto-throttle includes a feecjback of 
velocity error and of forward acceleration. The equations of motion of the 
auto-throttle system including engine lag are given by; 

'^th “ ^th 

SB 

e . = -0.05 - 0.5 a 

tn X 


where: 


(A. 13) 
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OF. POOR QUALITY 


For nanual throttle control is directly connected to the throttle. 


so that - 6^^. 


In order to compute the vehicle path, the components of V in the body 
coordinate system, have to be transformed in the components of V 

in the inertial coordinate system x., y., z. according to: 


y.‘ == • ^ 1 - 


where E. . = [F. * 


(A.IA) 


(Ael5) 


is the transformation mtrix for rotation from body to inertial system, which 
the transpose of the Euler matrix for rotation from inertial to body 

system. The Euler matrix is composed of a successive yaw, pitch and roll 
rotation, in that orde" according to: 




(A. 16) 
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(A.17a-c) 
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• • • 

The vehicle path is obtained by integration of a:*, y> and z.* 

V X- X 

The relation between the body rotations p, <7, r and the Euler angle 

• • • 

rates tp, 0, is given by: 
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(A. 18 ) 


and the inverse relations are given by: 

<p = p *f a tan 6 sin (p + r tan 9 cos tp 

0 *= q cos tp - r sin tp (A.lQa-c) 

ij; = r(cos tp/cos 6) + q(sin tp/cos 0) 


• • • 

The Euler angles are obtained by integration of tp, 0 and For 
relatively small values of 9 and tp, Eqs. (A. 19 a-c) become: 


tp = p 

0 = q (A. 20 a-c) 


= r 
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B. PPJIDICTOa LAGS 


B . . Prcd ctor Laws Based on a Circular Predicted Vehicle Path 


The predicted vehicle path is assuned to bo tangential to the velocity 
vector of the vehicle, I'. The projections of the predicted vehicle path on 
the x.ti. locally level plane and on the s.s. locally vertical plane, are 
hereafter referred to as the lateral and vertical vehicle path, respectively. 
The vehicle path is conjiuted frora the lateral and vertical path accelerations, 
vhich are derived hereafter. 

The translator)' equations of motion of the vehicle arc given by: 




1 .. n 
» — A * a 


i*. ^ Ki-, « r«. - »7 cos 0 sin — T ft 
b b b. ' n 

» 1 r? 

"’l. ♦ ~ cos 'P 

:> '? 0 n "h 


where irl , dy] and .il constitute the specific forces measured by nccclero- 

meters in and a, ^ body axis directicn, respectively, Tlie inertial 

accelerations in body axis direction are given by: 


a IK ♦ QW, . 

^ 




iron rqs.(A,ri) and Lqs . it follows tJuit the inertial accelerat i ons 

in hodv axis direction can be computed from the measured acce 1 erat ions by: 
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o' » a” ~ g sin © 
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^ cT ♦ ^ cos 0 cos U> 
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(A.23a-c‘) 


The lateral and vertical path accelerations and are then obtained by: 


i 1 • 

a, = 'i cos O - a, sin tp 

^ yb 


a - u sin cos 


(A.24) 


(A.25y 


Assuming the lateral and vertical path accelerations remain constant over 
the prediction span, a circular lateral and vertical vehicle path is obtained. 
TJie instantaneous lateral and vertical path-angle rates X ^ co.iputed 


XCt) = and = < 2 y(tyi' 


(A.26a,b) 


respectively. The situation for the horizontal path is shorn in Fig. 3 j. The 
actual path-angle and lateral deviation are denoted by X and respectively.and the 
predicted path-angle and predicted lateral deviation by Xp and respective!) . 

For relatively small angles of <>. B and X. the predicted path angle and 
lateral deviation, T seconds in the future, are given by; 


X.,(^) = X(t) ♦ 


(A. 27) 


tt) = !.•,;(:) + X(i) 


(A.2SJ 


1 + vmAv( 


vCr) = 2’ xCt) 


(A. 29) 
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and D is the predictor distance for a prediction T seconds in the future. For 
constant V the distance D is given by i? = TV, Note that the lateral 
displacement of the vehicle symbol on the display screen is given by the 
visual angle at which the symbol is seen in the body coordinate system. 

The lateral visual angle is given by: 

KCt) = xct) ♦ ^ (A.30) 

liqs. (A. 29) and (A.30) shows that the variables necessary for positioning the 
vehicle symbol on the display, are the path angle and path-angle rate only. 

For small y, actual path-angle and lateral deviation are related by: 


= xCt)/!' 


(A. 31) 


Substitution of Eqs. (A. 31) and (A. 29) in Eqs. (A. 27) and (A.2S) yields after 
1 inear lication: 


Xp(^) = ^ xCi) (A. 32) 

^ 

p 


which are Taylor series expansions, until the first and second derivative. 


Similarly, the predicted bank-angle can be approximated by: 


cp^Cr) = <p(t) + T (A. 34) 

The circular predicted vehicle path given by Eqs. (A. 27) and (A. 28) is 
useful in particular in following constantly curved sections of the trajectory. 



In the steady state, the radius of the vehicle path will be equal to the 
radius of the trajectory which enables a zero steady-state error in the 
lateral deviation. 

The equations for the vertical vehicle path are identical to Eqs. (A, 27) 
and (A. 28). Since constantly curved sections do not exist in the vertical 
profile of the trajectory, the linear approximations of Eqs. (A. 32) and (A. 33) 
will be sufficient for the vertical vehicle path. 


B.2 Predictor Laws for Forward Velocity Control 

For a given prediction time T, the predictor distance D can be computed 
from the present states of the longitudinal dy^namics and from the present 
control commands. Since the engine dynamics and velocity control are of a 
much lower natural frequency than the short-period d>Tiamics, the computation 
of D can be simplified. The engine dynamics are represented by a first-order 
lag with time-constant 1/c seconds, given by: 


2 = -ea + k 6. 

X X th 


(A. 35) 


where is the foiward acceleration, due to engine thrust and 5^^ is the 
throttle position setting. The fonvard velocity U is the integral of 
and the travelled distance S is the integral of U. The state equation for the 
simplified velocity control system is given by: 
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(A. 36) 
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or vnritten shortly as: 

= '■t ?t<« * h 

Assuming 6^^ remains constant over the prediction span, a linear 
T seconds in the future, from the present time i onwards is given by; 

x^(t+T) = <I>^(T) a^Ct) + r^Cx) 

where (x) is the transition matrix at T and 


where 


and 
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(A. 37) 
prediction, 

(A. 38) 
(A. 39) 

(A. 40) 


(A. 41) 


The travelled 'distance at (t+x) is given by: 





ORIGINAL IS 

OF POOR QUALITY 

- 32 - 


5(t+T) = ^ + ct l) a^it) + T Uit) + 

" + <?T - 1^ (A. 42) 

For T = T the predictor distance Z? is given by: 

D = - Sit) (A. 43) 

Expanding the exponential term in Eq. (A. 42) until the third derivative and 
substituting Eq. (A. 42) in Eq. (A. 43) yields: 





ajt) 



it) 


(A, 44) 


Tnc vehicle s>Trbcl with velocity tick-marks will be positioned at distance 
D ahead, whereas the solid reference square at distance Do - 'T Uo ahead, 
where J: is the desired velocity, see Fig. 1. 'fhe velocity U is brought to 
the desired value *‘r by bringing the distance between D and Do to zero. Thus 
the error for velocity control is defined as: 



(A. 45) 


and the structure of the control system is shown in Fig. 34a. 


The effect of variations in U are strongly observed in apparent changes in 
sire of the predictor 5 \mbol. Ill us , an increase in L yields an increase in D, 
which is manifested in an apparent shrinking in size of the vehicle svanbol . 

On the other hand, a decrease in U yields an apparent increase in size of the 
veh id t‘ s vnhol . Fq . 



(A. 44) also shews than an increase in engine thrust due 
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to a forh-ard throttle displacement, results in a simultaneous forviard notion 
of the vehicle symbol, so that throttle displacement and vehicle s>Tnbol notion 
are in the same direction. 

Note that, apart from a change in apparent s>Tnbol size, the location of 
the vehicle s>Tnbol on the image is not affected by variations in predictor 
distance. This location depends on the visual angle only which does not 
depend on Z?, see Eq. (A.30). 

Tiie diagram of Fig, 34a is further simplified in Fig. 34b, which shows 
the existence of a second order lead-term in the feedback path, given by: 


Dis) 





CA.46) 


Tlie zeros of Eq. (A. 46) have a natural frequency which is inverse proportional 
to T according to u' p 2,45/r rad/scc and a fixed damping ratio of Z = 0.61, 
Fir.cc the pilot response is much faster than the system response, the pilot 
transfer function can be considered as a simple gain K . For sufficiently 
large values of K^, the closed loop system poles will be located at the zeros 
of £q. (A. 46), and will be the dominant system poles. A small value of T 
yields a well damped, fast response, whereas a large value of T yields a 
sluggish response. However, a fast response will be at the account of 
considerably increased throttle activity. A value of T should be chosen which 
yields a satisfactory response with still acceptable throttle activity. 

Since the dynamics of the velocity control system are of a much lower 
natural frequency than the lateral and vertical d>Tiamics, the value of T, 
optimized for velocity control will probably be too large for the lateral and 
vertical coTitrol . Fig. 34c shows tbiC realization of a system in which the 
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prediction time for velocity control differs from the prediction time for 
lateral/vertical control. The prediction time for velocity control is 
denoted by and the prediction tine for lateral/vertical control is 
denoted by T, Like in the previous cases the solid reference square is again 
positioned at Do - T Uc. However, the deviation d = Do - D between the 
solid square and the tick-marks, is computed by means of prediction time T^, 
Adequate results are obtained for - 2T , 


B.5 Full-State Predictor Laws 

Tlie full state predictor law is derived for the lateral predictor notions 
only. The derivation is valid for the vertical predictor notions as well. 

The lateral dynamics Avith stability augmentation system engaged, are derived 
in Appendix A, I*or simplicity, the yaw -rate wash-out filter and ^ervo 
activator dynamics are replaced by unity gains. The augmented lateral system 
is given by: 
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(A. 47) 


where is the system matrix, B, the control input iratrix and U\ the 
disturbance input matrix of dimensions (r; ’^ »:) and , respectively, 

and 
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For small angular deviations from a straight reference trajectory, the 
hcadinv\ angle 'j-' and lateral deviation are obtained by integration of. 

* = (A. 49) 





: ? 

4 
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and = 7(\}i *► B) CA.50). 

Addi»'<T Eqs.(A.49) and (A. 50) to Eq. (A.l) yields the following augmented system: 


• 

s 


1 

1 

• 

p 


1 

■ 


Ao 1 

tp 


1 

• 

I* 


1 

i 

'I' 


T 

1 

0 0 0 1 1 

1 

k 


I' 0 0 0 1 


0 0 
V 0 




Bo 


Wo 


V 

-9 


(A. 51) 


Defining 

= col[i!;,^^] 

Eq. (A, 51) can be written as: 


(A. 52) 


• 

xo 


Aq 

I01 


Xo 


Bo 


Wo 


_ 

f 



+ 


u + 


• 

f 1 


Pi 



XI 


0 


0 

— 


k 1 -J 


— 


— 


— 


V 

-9 


(A. 55) 


Inhere Pj and P 2 are partition matrices of dimension (2 x m) and (2x2), 
respectively, Eq. (A. 53) can be written shortly as: 


ar = i4r + Pu + Wv 

- - - -9 


A 


ulierc r = 


LI J 


(A. 54) 
(A. 55) 


‘•5 
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A linear prediction T seconds in the future, from the present tine t 


onwards, is given by: 


ft t f'T 

x(t+T) = <I'(T)x(t) + ] 4-(T)B«(t-4T)d0 + J HtW‘>git*r)da (A.56) 


where 4>(t) is the transition matrix at x of the system of Eq. (A. 54) . The 
prediction of Eq. (A. 56) consists of three parts: the first tern at the 

right-hand side constitutes the homogeneous system response, the second 
term constitutes the forced motion due to control commands and the third tern 
constitutes the forced notion due to external disturbances. 

Assuming M(t)- and y^(t) remain constant over- the prediction span, 

Eq. (A.56) can be simplified as: 

x(t+T) = <t>(x)r(t) + r(T)Bw(t) + r(x)lv'u^(t) (A. 57 ; 

where r(x) is the first integral of 5(x) according to: 


'(T) = 


I «(c)cfa 


and r(T) are given by the series expansions: 


4(t) = J * At * Y 


a" 4 

nl 


o T ^ ' T 

r(T) = IT * A yy * ^ jT “ * • • • 


Expanding the system matrix A by means of Eqs. (A. 59) and (A. 60) yields: 
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0(T) = 


and 


rcT) = 


<5>o(T) 

0 

PiEoCx) + P2 PiAo(t) 1 

I + PzT 

ToCx) 

1 0 
! 

PiAo(t) + PaPiaoC*!^) 

1 

1 It + P2 


(A. 61) 


(A. 62) 


where is the transition matrix of Ao and To, Ao and «o are the first, second 
and third integrals of respectively and are given by the following series 
expansions: 


2 n 

^0 C"^) = J + AqT + Ao yT **' • • • ^ + • • • 


2 3 

To (T) = It + jr ^0 3T * • • * 


n~\ 

n\ 


+ . , 


T* , T** ,2 T® .«-3 t” 

5o(T) =Ijy- + i4ojj-+'4o jj- + ... + >5o + ••• 

» 

From Eqs. (A. 63) - (A. 66) the following relations can be derived: 
_ 

A(j(T) = Ao-o(T) + I 
ToC") = /1 oAo(t) •*• It 
’’cCt) = /iofoCt) + I 


(A. 63) 
(A. 64) 

(A. 65) 

(A. 66) 


CA.67) 
(A. 68) 
(A. 69) 


Substitution of Eq. (A. 68) in the lower left partition of Eq. (A. 61) 


yields: 



«(T) = 


Oo(T) 


PjAoCT) + PiT J + PzT 
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(A. 70) 


and substitution of Eq. CA.67) in the lower left partition of Eq. (A.62) yields: 


ToC-r) I 2. 

rcx) = 


PaEoCx) + Pi - 5 - , IT + Pz 


CA.71) 


where P 3 is a matrix or 


' dimension (2 >|m) given by: 


P 3 = PiAo P 2 P 1 


(A. 72) 


Evaluating Eq. (A.57) with Eqs. (A.55) and Eqs. (A.70) and (A.71) yields the 
prediction of ri : 

Xj(t+T) = P3Ao(T)Xo(t) + PlTXo(t) + (X + P2T)X1 (t) + 


+ P 3 Ho(T)Bow(t) + Pi ^Bou(t) + 
♦ p^^c(J)»oVg(.t) + Pi — 


(A. 73) 


Provided xcC^) and xi(t) are known exactly, the first terms on the right- 
hand side of Eq. (A. 56) can be computed exactly. However, the remaining 
terms depend on the future values of u and which are not known. Kith 
respect to the future values ''f u and v^, three assumptions can be made which 
arc shown in Fig. 35: 

(1) u and renain constant over the prediction span; (2) u and 
are set to zero; (3) u and v decay exponentially over the prediction span, 


X 
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according to: 

u(t+T) = M(t) 

u (t+T) = V Ct) 

0 0 
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(A.74) 
(A. 75) 


where • and are the time-constants of the decay. It is clear that x^^ 

or x^ identical with assumption (1) and also that x^^ or x^^ •♦■0 is 

identical with assumption (2). Assuming u(f) and V C^) can be approximated 

if 

by first-order Markov processes with spectral densities: 


4> 

uu 


(SO 


1 

(fli-s) (ai*s) 


*^u V ~ (G2-S) (<32+S) 

ff 9 


(A. 76,77) 


then the optimal prediction according to Wiener-Hopf is given exactly by 
Eqs.(A.74) and (A. 75), where x^^ = l/aj and 1/^2 . Thus the time-constant 

of decay is determined by the break frequency of the Markov process. 
Substituting Eqs.(A.74) and (A. 75) into Eq. (A. 56) yields: 


x(t+x) = «(x)x(t) + r*(x)Bu(t) + r**(x)f.^^(t) (A.7S) 

where F* and F** are expanded according to Eq. (A. 60) but with the matrix A 
replaced by 

A* = [A - aji] (A. 79) 

and A** = [A - C 2 JI (A. 80) 


respective!'' . 

The prediction of n, subject to decaying control and disturbance inputs, 


is then given by: 
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Xl 


(t+T) = PsAoWXot*) + PlTXo(t) + (I+PzT)Xi(t) + 

^2 


(A. 81) 


3 “0 0~g ' 

where the and superscripts indicate that the corresponding -xtrix is 

generated with: 

(A. 82) 


ard 


Ao = [4o - <21-^1 
ylo* = [Ac - C2I] 


(A. S3) 


respectively. 

1„ Appendix B.l it is sho«n, thet linearization of the circlar path viclds 
3 Taylor series expanslo". “'"•‘1 * prediction of ;i based 

a Taylor series expansion is given by: 


on 


X, (c) = Xl(t) -p TXl(t) + -^flCt) 


(A. 84) 


Differentiating the lower half of Eq. (A. 53) with respect to r and substituting 
in Eq. (A. 84) and using Eq. (A. 72) yields: 

Xl (t) = (r+TP 2 )XlCt) + PlTXo(t) + yrP 3 X 0 (t) + 

- ‘P ^ *“ * 

. §ir,Siu(e) 

The difference betueen the two prediction methods is defined as; 

A a y h ' (A.S6) 

ASi(t+T) = X-,{t^X) - Xi^(t) 

Subtracting Eq. (A. 85) from Eq. (A.Sl) and substituting Eq. (A.67) yields: 
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AriCt*T) = r,/1o=o(T)re(t) ♦ FJHJ (T)So^{t) ♦ 

♦ (A.S7) 

Eq. (A. 87) yields predictions for the heading angle and the lateral 

deviation, y ,, However, it was found wore useful to provide the lateral 
(1 

path-angle X rather than the heading angle, 5«o that the vehicle sn't^oI Is at 
all tines tangential to the vehicle path* A combined linear/non-lincar 
predictor law can be formulated by replacing the term in Eq. (A.S4) 

by a circular path. At t 7 the predictions for path-angle and lateral 

deviation at'c given by: 

xCt»r) “ x(i) * ?x(f) ♦ (A.ss) 


. I* 


rx(t) 




(A. 89) 


Eqs. ifV.88) and iA.S9) are identical with liqs . (A. 27) and (A.2S) of the circular 

path prediction with the addition of the terns :.x{t*T) and The 

» ' 

.i.idition for the lateral deviation follows iron Eq. (A. 87) and is given b> . 




(A. 90) 


wliore 1 . • • ] der.olos tlio second row of matrix {...). Hie addition tor the 
path-angle i- ot'tained In* using the relation: 
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(A,91) 


Differentiating Eq. (A. 87) with respect to T and substituting T = 7 yieHs: 
AxCt+r) » 1 jrowalPj/loAoCTlxoCt) + 

+ rowjPaA^CDBolwCt) ♦ 

+ row2[Pj’‘Aj'‘(3’)K'Jy^(t)| (A.92) 


The linear/non-linear predictor law given by Eqs. (A.SS), (A. 89) and 
Uqs. (A. 90) and (A. 92) combines the advantage of a circular path predictor law, 
i.e., the zero steady-state error in constantly curved sections, with the 
advantage of the sophisticated full-state predictor law of providing a more 
accuractc prediction. 

Tlic prediction for the bank-angle follows directly from Eqs . (A./S), (A.SS) 
and (A.48a) and is given by: 


- roWj[-.*'/r)]x5(t) ♦ row^ (r^''(r).Vj]r^(t) 




(A. 95) 





-0,0191 

O.OSIO 

0.0906 

0.0360 

0,0049 

-0.0134 

-0.0254 

-0.0095 

-0.1766 

0.4711 

O.S5.';9 

0.3324 

-0.0224 

0.0507 

0.1059 

0.0421 
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0.6851 

0.2870 

0.5523 

o.oooT 


-0.1766 

0.4711 

-0.1641 

0.5324 

To = 

-1.2420 

2.7950 

5.4960 

1.9040 


0.8008 

0.2864 

0,5871 

0.2663. 


“ 4.027 

0.772 

1.494 

-0.495 


-1,242 

2.795 

-0.504 

1.904 

It 

O 

C 

-3,815 

7.930 

17.01 

5. 187 


4.441 

0.653 

1.564 

1.031 





„ -T 


”ll.30 

1 .577 

2.775 

-1.723 

Co = 

-3.815 

7.930 

-0.9S7 

5. 1ST 


-7.466 

14.91 

54. SS 

9.544 


_ll.S3 

0.968 

2.76S 

2.779 ^ 

Tbo matrix is computed with l:q, 

. (A. 72) 

and is gi\ 


3 , 869 

-0.177 

0.0 

-3.377 

?3 = 






j^73.46 

0.0 

32.19 

46.36 ^ 

hl)cn assumim: 'S 

' CL 

remains 

constant over 

the prediction spe 

- 0 so that >»o = arid substitutini: the previous val 

and r.q. (A. 91) y 

Lelds : 




A;/ ( t + (> ) - 

1110. S + 

228. S p - 69.0 

ip - 583. 4 

r - 337 . 0 

and 





A>;(r >6) =1 

..SO B ♦ 0. 

337 r - 0.115 o -0.S40 

»' - 0.737 V 


(A. 9-:) 




The tivst order ro 1 1 -prediction of Tq, CA,3*)) is piven by: 
cp (.r) = 1 .0 (P 6.0 r 


(A,9?0 
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and the full-state roll prediction of Eq. (A. 93) is given by 




<p(t +2') = -0.1766 3 -f 0.4711 p + 0.8359 o + 0.332 r - 2.136 6^ (A. 97) 

Fig. 36 shows the coefficient of 6^ in the conputation of £>y^, as a 
function of T and for various values of Fig. 36 shows, that the curves 

can be closely approximated by quadratic functions, which might considerably 
simplify on-line on-board computations. Fig. 37 show's these cur\’es for the 
coefficient of 6^ in the computation of X- Also the curves in Fig. 37 can be 
approximated by linear or quadratic functions. 


C. ALGORITHMS FOR COMPUTING THE TRANSITION AND INTEGRAL RATRICES 

The computation of the transition matrix and integral matrices starts with 
the computation of the third integral matrix Ho(t), by means of the series 
expansion of Oq. (A. 66). Tlie expansion is continued until the norm of the n-th 
term is sufficiently small. The second integral matrix Ao(t) is obtained by 
substituting ao(^) into Eq. (A. 67), the first integral matrix To ( t) is obtained 
by substituting (t) into Eq. (A. 68) and finally the transition matrix ^oC'^) 
is obtained by substituting TcCt) into Eq. (A. 69). 

Computational difficulties may arise for large values of t, resulting in 
overflow before the series converges. For the computation of vs(i) this 
difficulty can be overcome by subdividing t in /•' intervals At according to: 

T = .VAi, and using the relation: 

45(T) S tc(.VAt) = {v\(At)}^‘’ 


(A.9S) 
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However, this relation is not valid for the integral natrices, and the correct 
expressions are derived hereafter. 

Given the homogeneous system equation of order n: 


2 q = AzZ^ 


The first, second and third integrals of 2c are defined as: 


‘ A • A . • ^ 

:i = 2o J Z2 - , 23 - j-2 


and an augmented system of order is defined as. 


21 


23 


i i- 


^ „ 1 / 


23 

21 

i 5-’ 


(A. 99) 


(A.lOOa-c) 


CA.lOl) 


or shortly written as: 


(A. 102) 


Hie transition matrix is obtained by expanding A according to Eq. (A. 59), 
and is given by: 


rt^(t) 10 1 0 I 0 j 

j L-J- 4— 4— ; 

i r-(-) I I I 9 j 9 i 

” i -.(') i Ji 1 J i 0 I 

1 I ^ ^ \ 

, 1 . j 


(A. 103) 


= E;(T) 


1 


1 
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where the first partition column contains the required matrices. The relation, 
given in Eq. (A.98) is applied to the augmented system matrix of Eq. (A. 103) 
and thus $(At) is raised to the power N by (//-I) repeated multiplications. The 
u-th multiplication yields a matrix, from which the first column partition 
contains the following expressions: 




= r”'5'n + r. 

3 0 0 0 


♦ A 

G 0 0 0 0 


-C -0'*'0 ■*'9! ^G “0 


where: 


0 " 0 ’ 


t«Ai) ; r” r.OiAx) etc. 


0 0 


and 


= ; r, etc 


0 0 


(A. 104) 
CA.105) 

(A. 106) 

(A. 107) 

(A.lOSa.b) 

(A.109a,b) 


The algorithm starts with subdividing x into A’ sufficiently small intervals 
A" . Next, 5o, Ac, Fc and 4>ci are computed by means of the expansion of Eq. (A. 66) 
and by means of F.qs . (A. 67), (A.6S) and (A. 69). These matrices are substituted 
in Eqs. (A. 104) - (A, 107) to obtain the value of the matrices at n = 2 or (2Ax) . 

Tlie results are used in the next iteration to obtain the values at m = 5 
or (5ix) . Thus, starting with »: = 1 the process is repeated until n = .V- 
which finally yields the values of the matrices at x = .Vix 


1 
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D. ca 4 PirrER-GRAPHICS METHOD FOR DRAWING THE 3-D PERSPECTIVE VEHICLE SYMBOL 

A vehicle symbol, which is generated as a -’wire-frame" structure is 
ambiguous because it is not clear x,rhether the vehicle symbol is pointing 
towards or away from the observer. In order to obtain an unambiguous vehicle 
symbol, the hidden lines should be removed from the image. However, commonly 
used algorithms for removing hidden lines, require considerable computational 

efforts which might unfavorably affect the update rate. 

An alternative fast method has been developed to remove hidden lines from 
the image without requiring additional CPU time. Although the vehicle s>'mbol 
chosen in this program has a square and constant fuselage cross-section this 
method applies to any. not necessarily constant cross -section. The hidden- 
line removal is restricted to the fuselage structure. Thus the fuselage 
appears "solid" but the wings and tail still appear transparent. In spxte of 
these restrictions a clear and unambiguous image is obtained with minimal 
computational efforts. 

me basic vehicle symbol structure is defined by a set of coordinates in 
the object coordinate system between which straight lines are 

drawn, see Fig. 38. The coordinates are given in Table VII. The image is 
observed in the eye coordinate system (aligned with the body coordinate system) 
with the observer's eye at its origin, see Fig. 39. Thus each coordinate 
in the object coordinate system is transformed into a coordinate in the eye 
coordinate system according to; 




— 




or 1 
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e 


(A. 110) 
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where dx^, dy^, dz^ is the location of the center of gravity of the vehicle 
symbol in the eye coordinate system and 

Ve " (A. Ill) 

^ P P P 

is the Euler matrijf for rotation from object to eye coordinate system given 
by Eqs . (A. 16) and CA.17a-c), The sjTnbol image is obtained by projecting the 

coordinates on an image plane in front of the observer and perpendicular 
to the axis. The limits of this image plane are from +45 to -45 degrees 
horitontally and vertically. In order to draw the image, the coordinates on 
the image plane have to be interconnected by lines. The pattern in which the 
coordinates are connected is defined by a list of move/draw instructions. 

Each instruction is a coordinate number, which can have a positive or a negative 
sign. A positive sign indicates that the CRT beam moves from its previous 
location to the specified coordinate whi le drawing a visible line (draw) and a 
negative sign indicates a shift to tl\e specified coordinate with the CRT bean 
blanked (move). Ilic first location in the list contains a count of the total 
number of move/draws. Special treatment is required for lines from which one 
oi both coordinates fall outside the limits of the image plane. In that case, 
the line is only partially visible and the intersections of these lines with the 
image plane borders have to be computed (clipping). An efficient algorithm 
is implemented for handling this problem. 

In order to determine which lines of the predictor fuselage arc visible, 
the position of the observer's eye in the object coordinate system is to be 
computed according to: 
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dx 




dx 

P 




e 

% 

= - 




dz 





P 




e 

L- - 

- 

*— - 




(A. 112) 


where E is the Euler matrix for rotation from eye to object coordinate 
system and 


E = {e y 


(A. 113) 


e'^ p-^ 

This transformation is shown in Fig. 39. The minus sign indicates that if the 
displacement of the vehicle position in the eye coordinate system is positive, 
the displacement of the eye position in the object system should be negative. 
Since E was already computed previously, its transpose is readily availabe 

p-tB 

and the wansformation of Eq. (A. 112) can be carried out with minimal computa- 
tional effort. Seen from behind, the observer's eye can be located either 
above or below each one of the four planes applied through the fuselage sides, 
see Fig. 40. These planes are given by: 

= b (A. 114) 

= -b (A. 115) 

= b CA.116) 

= -b '1.117) 

where b is half the diagonal of the fuselage cross-section. Altogether, there 
are nine areas in which the eye position can be located. In each of these 
areas different lines of the fuselage will be visible. A four-bit code 
determines the area in which the eye position is located. In order to compute 


for 

right upper side 

A : 

Z + V 

p p 

for 

left lower side 

C : 

p p 

for 

left upper side 

D : 

z - u 
p '■p 

for 

right lower side 

B : 

- ' 

p p 
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this code, four tests are made and accordingly, the appropriate bit in the code 
is set in the following way: 

above righ-upper : dz^ + dy^ > b bit 0001 is set 

below left-lower i: dz + dy < -b bit 0010 is set 

\ P P 

above left-upper : dz^ - dy^ > b bit 0100 is set 

below right-lower: dz^ - dy^ < -b bit 1000 is ^et 

For each one of the nine areas, an individual inoYe/draw list is defined. 
wnich only draws the lines visible in that particular area. All the individual 
move/draw lists are stored sequentially in one area of the memory. 

According to the value of the code, the appropriate move/draw list is 
referenced through an address table which, for each code, contains the address 
of the first location of the required list. If the code remains cero, none of 

the fuselage sides are visible and the image is completed after drawing the 

always visible parts of the image, e.i., wings, tail and rear surface of 
fuselage. A flow-chart of the computational method is shown in Fig. 41. 

Tne movc/draw lists are given in Table VIII, and the address table in Table IX, 

The method can be applied equally well to a fuselage with a non-square 

cross-scction. For example, for a triangular cross-section the total number 
of areas, in which the eye position can be located is 7 and each area is labelled 
bv a 5-bit code. For a hexagonal cross-scction thcrearc 19 possible areas, 
which are labelled by a 6 bit code and for an octagonal cross-section, there are 

55 possible areas, labelled by an 8 bit code. 

The vehicle s^bol with hexagonal cross-scction, is defined in Fig. and 
the coordinates given in Table X. Ihe possible areas in which the eye position 
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can be located, are shown in Fig. 43, In order to compute the 6-bit code, six 
tests are made and the appropriate bit is set as follows: 


for 

dz + — dy 

P P 

> b 

(above A) 

bit 

000 

001 

is 

set 

for 

dz + — dy 

P P 

< -b 

(below D) 

bit 

000 

010 

is 

set 

for 

dz — t dy 

P Vz P 

> b 

(above F) 

bit 

000 

100 

is 

set 

for 

dz — t dy 

P Vz P 

< .b 

(below C) 

bit 

001 

000 

is 

set 

for 

dy > Ir'Jz b 
'P 2 

(right of B) 

bit 

010 

000 

is 

set 

for 


(left of E) 

bit 

100 

000 

is 

set 


where b is half the diagonal of the hexagon. 

The codes for the 17 areas are shown in Fig. 45, and the corresponding 
move/draw lists and address table are given in Tables XI and XII respectively. 
It is clear that a more complex fuselage cross-section will primarily increase 
the storage requirements for the move/draw lists, and, to a much lesser extend 


affect CPU requirements. 
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TABU: I. Rr.SULIS OF IRAJLCTOkY FOLLOWiNd: SUIUIXT A. 
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"Strairht” Tunnel; Width: 300 ft; Foil Version 1 

"Banked" Tunnel 

2-0 Predictor 

3-0 Predictor | 

Manual Throttle 

3-D Predictor 
• ISOO ft 

fj • 
300 ft 

p,‘ 

ISOU ft 

Dq* 

2000 ft 

^0* 

900 ft 

D,‘ 

1500 ft 

Oo» 

2000 ft 

npull*' 

.State 

Predictor 

Uiw 

A-1500 ft 

3-n Predictor 

"tick" 

Marks 

U IV 

Digital 

Width: 
300 ft 

Width; 
450 ft 

Width: 
300 ft 

Read-Out 

Roll-Ve 

rslon 

Roll 

Stabilized 

NUMur.K ni 
RUNS 

6 

6 

8 

6 

6 

8 

6 

6 

b 

8 

6 

6 

Cov (n ; ) 

o.in 

.♦0.045 

1.008 

10.122 

1.943 

10.241 

0.250 

10.036 

1.213 

10.093 

2.053 

♦0.246 

0.972 

10.118 

0.857 

±0.216 

0.789 

±0.077 

0.623 

±0.060 

0.484 

10.064 

- 

Cov ) 

lin"ft‘V 

0,312 

♦0.0(»O 

2.001 

10.278 

4.066 

10,649 

0.759 

♦0.067 

2.742 

±0.187 

S.IOS 

10.63S 

2.715 

♦0.391 

4.269 

±1.737 

3.064 

±0.591 

2.010 

±0.199 

2,137 

±0,261 

- 

Cov CO) 

IioTt') 

0 JOO 
10 V 069 

0.361 

♦0.046 

0.438 

10.073 

0.316 

10.065 

0,288 

♦0.058 

0.488 

±0.081 

0,381 

±0.031 

0,318 

±0.048 

0.458 

±0.159 

0.374 

±0.057 

0.342 

±0.063 

0.SG3 

±0.111 

Cov (3,) 
lio’ft') 

0.331 

♦0.072 

0.561 

.♦0.123 

1.093 

10.364 

0.677 

±0.184 

1.042 • 
10.192 

1.803 

♦0.507 

1.270 

±0.367 

1.525 

±0.348 

1.198 

±0.259 

1.739 

±0,557 

3.461 

±0.802 

0.926 

±0.230 

Cov (1) 

fradM 

0.H2 

10.003 

0.134 

iO.OOl 

0.123 
±0.003 j 

111101011 

0.136 

±0.003 

0.123 

±0.003 

0.134 

10.003 

0.137 

±0.003 

0.141 

±0.003 

0.129 

±0,004 

0.126 

±0.003 

0.139 

±0.003 

Cov (f) 

(1 O' '(rad/ 
^ cc 1 ^ 1 

0.814 
? 0.037 

0.439 

♦0.056 

1 

0.468 

♦0.175 

0.565 

±0.105 

0,410 

±0.082 

0.368 

±0.053 

0.463 

♦0.062 

0.442 

±0.028 

0.456 

10.050 

0.386 

±0.063 

0.333 

±0,035 

0.534 

±0.077 

Cov (u) 

lio’d't/ 
•■•cc ) M 

- 

- 

- 

- 

- 

- 

- 

3.860 

±5.033 

2.892 

±1.184 

- 


- 

iiriJ i 

Cov 

110 (ft/ 
sccMM 

- 

- 

- 

- 

- 

- 

- 

0.686 
j ±0.096 

1.285 

±0.918 


- 

- 

Cov (*Vi) 
I10‘*rad^l 

1.102 

10.116 

0.901 

10.178 

0.642 

±0.334 

0.902 

10.216 

0.887 

±0.222 

0.450 

±0.112 

1.088 

±0.148 

1 0.582 

±0.077 

0.439 

±0.107 

0.444 

±0.099 

0.520 

±0.141 

1.120 

±0.152 

Cov (Ar) 

110''r.id'| 

0.332 

jn.oBfl 

0.216 

•0.051 

0.280 

10.136 

0.246 

♦0.060 

0.209 

♦0.035 

0.192 

10.045 

0.441 

±0,071 

0.310 

±0.118 

0.218 

±0,083 

0.199 

±0.035 

0.253 

±0.068 

0.237 

±0'.060 


ORIGinAL F/tGii iS 
OF POOR QUALITY 



tabl; n. a TfLurcioRr jcu.oi^inc: R'BJLCt e. 


?*t' rrcJictor 


"StrjiitS!" Turnrl; width: 100 ft, P>n Version 


,1-0 rrrdictor 


I o» 

1 l-'!AS 


'('JV (r;/. : 

Mi'j'ft' 1 

\ 

ll‘v)V I * l) 

illM'l! ' I 

!fov 

I 

I lln’ft* ) 

( ; ) 

Mr»d’ } 

It MV tf ) 



\C’.U \^<} 

|{ln'(ft/ 



\i\>v f\;,) 

I U> i ft / 

jCJv 1) 

(lo'^radM 


V>i • 

ir»oo ft 

In • 

:coo ft 

! . i ‘ 

1 i.o:i 

2.1U 

* 

1 

*0.233 

1 : . t.v * 

S.4f>9 

j 

•0.824 


COO ft 900 ft ir»00 ft 2CH)0 ft 


1..V0 I o.:9'> 


iUnuJl Throttle 
3-D rr?4lctar 
Tv . 15C0 ft 


-lyj!” 

5l#tc 

Predictor ,, 

law , 

^I•1SOO ft 


•'tlcl'* DlslUl 

)*ark» Read-Out 



1.25: 

Ui.US 


^'Banked** Tunnel 


3*0 Predictor 
Dt • ISOO ft 


Rol 1 -Version 


0.670 
to. 061 

0.S96 

10.093 

- 

2.073 

10.2S4 

1,746 

10.22S 

- 

0.236 

10.028 

0.374 

10.062 

0.272 
to. 014 

0.974 

iO.MS 

1.010 
1 so. 192 

0.19& 

S0.1S7 

0.134 

10.002 

0.129 

10.003 

o.m 

10.003 

0.374 

10.027 

0.281 

10.023 

0.4S7 

10.064 


• 

- 

- 

- 

- ’ 

0.687 

10.142 

0.780 

10,160 

0.991 

10.144 

0.387 

10.133 

0.498 
10. OSS 

0.S39 

10.129 


OF POOR QUALH V 






TABU. in. Pr.SULTS OF TPAJrCTDPy rOLLCHiNC: scwca c. 



'*Straj;:ht” Tunnel; Width: 500 ft; Roll VcMlon 

*'BarJted” Tunnel 



i-D Predictor 

j 5-D Predictor 

1 Manual 1 

rhrottle 
‘dictor 
500 ft 

J-D Predictor 
D, . ISOO ft 

Di . 
900 ft 

T'8 . 

ISOO ft 

Oo » 

2000 ft 

Di ■ 

900 ft 

/;; • 
1500 ft 

t/'j • 

2000 ft 

*Tuir’ 

State 

Predictor 

Law 

/7,-1500 ft 

1 3-D Pr< 

i 

•Tick" 

Marks 

t 

OlKltal 

Retd-Out 

Width: 
300 ft 

Kldth; 
450 ft 

Width: 
300 ft 

Roll -Vers Ion 

Roll 

Stabilized 

MWBER or 
kONS 

6 

6 

6 

: 

6 

7 

6 

1 ^ 

7 

6 

6 

7 

Cov (i/^J 
(lo'ft'i 

0.0«5 

.*0.012 

0.568 

iO.087 

1.646 

10.229 

0.107 

20.048 

0.681 

20.048 

1.842 

10.249 

0.941 

20.113 

1 0.831 

1 10.057 

1.06S 

20,135 

0.652 

10.035 

0.716 

10.070 

- 

Cov ) 

0.430 

i0.042 

1.520 

10,327 

4.294 

10.812 

0.658 

J0.17S 

2.024 

10.285 

5.03B 

10.537 

2.396 

20.421 

1 3.543 

10.552 

3.027 

20.388 

2.563 

10.137 

3.368 

10.537 

• 

Cov 1 

(lO'ft’l 

0.313 

t0.07(> 

0.231 

10.027 

0.412 

10.154 

0.223 

10.046 

0.280 

10,068 

0.505 

10.120 

0.397 

10.048 

j 0.265 
10.060 

0.472 

20.124 

0.326 

20.046 

0.248 

10.024 

0.486 

10.072 

Cov (z^,) 

1.190 

m),h: 

0.969 

!0.151 

1.748 

10.611 

1.574 

20.391 

1.448 

20,278 

|ffil 


1.630 

20.359 

1.783 

10.539 

2.637 

10.192 

4.523 

10,783 

m 

Cov (^) 
[raJ'l 

0.153 

M).005 

0.136 

lO.OOS 

0.133 
! 10.004 

0.148 

20.007 

0.137 

20.002 

0.130 

10.005 


0.141 

20.004 


HQnHjli 

0.13S 

10.002 

0.133 

10.004 

Cov (D 
llO“?(rad/ 
sec)M 

1.018 

i0.12j 

0.746 

10.111 

0.743 I 
10.099 j 

0.769 

20.060 

0.S74 

20.057 

0.438 

10.059 

0.596 

20.092 

0.538 
1 20.053 

m 

0.402 

2C.061 

0.344 

10.044 


Cov fu) 

|I0'(ft/ 
ICC)' ) 

- 

- 

■ I 

• 

- 

- 

n 

1.390 

20.429 


- 



Cov C‘th) 
(lOCft/ 
sccM* 1 

- 

- 

- 1 

1 



- 

- 

- 

0.476 

20.123 

mmm 

- 

- 

- 

Cov (ia) 

[ 10“ * rad* ) 

0.962 

M).145 

0.704 

10.102 

m 



0.316 

10,097 

1.172 

10.264 

0.589 

10.119 

1.132 

10.297 

0.452 

20.087 

0.410 

10.072 

0.824 

10.141 

Cov (''t’) 
(lO-^rad^ 

o.2r.s 

*-0.078 

0.214 

10.031 

0.179 ■ 
10.035 

0,262 

20,043 

0,227 

20.077 

0.241 

20,061 

m 

0.262 

20.063 

0.319 

10.076 

0.199 

10.038 

0.206 

10.030 

0.187 

10.037 
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TAhi r !v, p.i.sui/iS / roK'.c’A’N*;: 




• *> 



i 


••BanVcd” Tunnel 


Manual Throttle 1 predictor 

5-D Predictor | Va • ISOO ft 

Do • 1500 ft j 

i ■ 

! 

j 

t 

1 

> 

i ■ 

> . ' 

'i 

1 1 : 

Ul 3 

1 M ■ ■ 

> ■ 

1 ; 

\ f 

\ 

i ! 

o o • 

j ,'j 

•0 c> f i 

8l^i 

:o l- . 

^■-r; i' • 

c: 7^ ‘ ■ ; 

j.' ; . ’.‘i 

C ^ 

^ (2 i 

' i J 

J 1 ■ . 

■•Tick" 

Marks 

Difiltal 

Read-Out 

Width: 
300 ft 

Width: 

450 ft i 

Bl 

Roll-Version 

Roll 

Stablliicd 

8 

6 

7 

- 

- 

0.771 

10,148 

1.161 

10.064 

0.670 

10.134 

- 

- 

1.241 

iO.251 

4.397 

iO.449 

2.650 

40.453 

- 

- 

0.301 
10. OSS 

0.470 

10.090 

0.345 

10.089 

- 

- 

1.574 

iO.472 

2.341 

10.371 

2.173 

10.336 

- 

- 

0.135 

10.005 

0.140 

10.004 

0.136 

10.004 

- 

- 

0.416 

10.072 

0.468 

10.078 

0.440 

10.097 

- 

- 

0.922 

0.243 

2.501 
1 .479 

- 

• 


0.851 

0.177 

2.447 

0.595 

- 

- 

- 

0.602 

iO.lOS 

0.899 

10.304 

0.546 

10.119 

- 

- 

0.327 

10.090 

0.271 

10.073 

0.237 

10.022 

- 

- 


< ' ; 

1 

< 


1 

I 


I 




































































”Straip,ht’' Tunnel, Width: 300 ft; Auto-throttle 


SMBJI'CT 


IMthlUCTOR 


No. of Scries 
of () Rims 


2 - 1 ) 


6 


j , 

;r jf 

[ 1 ()■' see 1 

r 

1 lO^’sec 1 

(io'ff-j 

r.c r.-:J 
(I 

Cov ('{’) 

jratlM 

I 

(lov ([ ) 

I ( Kf' (r;ut/ 

ii-mi 

Cov 

lio'* rad" J 

Cov (3c) 

1 lo’' rad" ) 


0.19‘J 

H).()14 

0,127 

10.006 

0 . 201 ) 

10.038 

0.249 
JO. (no 

0.197 

.♦ 0.010 

0.402 

10.080 

0.313 

10.06S 


3-1) 


6 


0.205 

n).027 


0.140 

♦0.014 


0.211 

n).054 


0.394 

♦0.040 


0.199 

10.010 


2 - 1 ) 


0 . 207 
♦0,009 


0,129 

10.004 


0.213 

n).019 


0.296 

10.02S 


0.197 

10.004 


0.365 

♦0.107 


0.263 

10.084 


0.331 

10.018 


0.226 

10.013 




3 - 1 ) 


0 . 206 
10.006 

0.142 


10.013 


0.216 

10.011 

0.56S 

10,022 

0.195 

10.002 

0.238 

♦ 0.020 

0.164 

10.013 


0.620 

10.0,53 


0,697 

10.174 


0.644 

10.053 


0.638 

10.049 


Roll-Version; Po - 1500 ft 


2 - 1 ) 


3-D 


2-n 


0.191 

♦0.009 


0,189 

10.010 


0.195 

lO.OlO 


0.125 

10.002 


0.132 

iO.Oll 


0.124 

i0.003 


0.163 

i 0 .v )18 


0.155 
to. 005 


0.175 

to , 016 


0.278 

♦0,023 

0,222 

iO.009 


0.391 

♦ 0.022 


0.225 
to. 003 


0 . 299 
♦0.037 


0.221 

i0.006 


0-. 494 
♦D.099 


0.406 

to . 022 


0.465 

to . 080 


0.330 
to. 064 


0.259 

10.016 


0.328 
to. 062 


0.659 
to. 04 2 


0,765 

10.076 


0.664 

10.053 


3-0 


6 


0.182 

±0.006 

0.125 

± 0.002 

0.167 

±0.011 

0.358 

±0.042 

0.214 

±0.006 

0.367 

±0.039 

0.244 

±0.027 

0.654 

±0.043 


I 


5 

\ 

i 



U1 

oo 


o o 

•n 




•c f > 
O ' 

a ^ 

TJ r'* 

O - 7 ) 
G V 


r- M 

3 k 
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OF FCGR Q'-'ALlTr 


TABLE VI. STABILITY DERIVATIVES OF DC-8 AIRCRAFT IN APPROACH-TO-LANDING 


Lateral 


= -0.113 sec“‘ 

= 0 . 

= 0.0238 sec Vrad 
= “1.328 sec’^ 

= “0.951 sec"^ 

= 0.609 sec ^ 

= “0.726 sec ^/rad 
= 0.1815 sec ^/rad 

- 0.757 sec ^ 

= “0.124 sec“^ 

= “0.265 sec"^ 

= “0.0552 sec'^/rad 
= -0.589 sec*^/rad 
= 243.5 ft/sec 

= 32.2 ft/sec" 

= 3,090,000 slug ft 

= 5,550,000 slug ft^ 

28,000 slug ft^ 


Longitudinal 

X 

u 

s; 

-0.0291 sec"^ 


c 

0.0629 sec"^ 


s= 

0. 

e 



z 

u 

= 

-0.2506 sec"‘ 


= 

-0,6277 sec“^ 

h 


“10.19 ft sec*^/rad 

e 



M 

u 

= 

0. (sec ft)~* 

u 

= 

-0.00S7 (sec ft)“* 

lU 

w 


i 

-0.0011 £t‘^ ! 

M 

= 

-0.7924 sec"^ 

1/ 


-1.35 sec"^/rad 

Uo 

= 

243.5 ft/sec 

Wo 

= 

0. 

9 

:= 

32.2 ft/sec^ 

^flap 

I 


35 deg. 






TABU: VII. UST OF COOUDINATHS OF PHRSPECTIVE VFHICLH SYMBOL Will! SQUARE FUSELAGE CROSS-SECTION 


, .. . 

Coorcji - 
nate 
Nu!r.!>cr 

Value in units of haJ f 
the fusclnjje width: d 

Desenp- 
ti on 

Coordi - 
natc 
Number 

Value in units of half 
the fuselage width: d 

Descrip- 

tion 

— 

r 

P 

V 

-p 

% 

X 

p 


z 

p 

1 

8.0 

1.0 

0. 


12 

-8.0 

2.S 

4.0 

hori- 










zontal 

2 

8.0 

0. 

1.0 


13 

-8.0 

-2.5 

4.0 

stabi- 

5 

8.0 

-1.0 

0. 


14 

-5.5 

0. 

4.0 

liter 

4 

8.0 

0- 

-KQ 

fuselage 

15 

4.0 

1.0 

0. 


5 

-8.0 

1.0 

0. 


16 

-3.0 

6.0 

0. 

left 










wing 

6 

-8.0 

0. 

1 .0 


17 

-3.0 

1.0 

0. 


7 

-8.0 

-1.0 

0. 












18 

-3.0 

1.0 

0. 


8 ■ 

-S.O 

0. 

I 

b 


19 

-3.0 

-6.0 

0. 












9 

-8.0 

0, 

4 .0 

vertical 

20 

4.0 

-1.0 

0. 

■m 

10 

-7.0 

0. 

4 .0 

.stabi- 










1 i zer 

-21 

0. 

0. 

6.0 


11 

-4.0 

0. 

1 .0 












1 22 

0. 

0. 

8.0 

vertical 






)} 




tick- 






23 

0. 

0. 

-6.0 

marks 






1 24 

0. 

0. 

-8.0 



{. 
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TABLE VIII. MOVE/DRAW CODING LIST FOR PERSPECTIVE VEHICLE S^TiEOL KITH 
SQUARE FUSELAGE CROSS-SECTION. 


Code 

(Binary) 

Code 

(Dec- 

imal) 

List 
Address 
in Core 

M/D 

Count 

Coding List 

Description 

0000 

0 

A 

22 

-6 , 7, 8, 5, 6, 

Rear fuselage surface 





9,10,11, 

vertical stabilizer 





-14,12,13,14, 

horizontal stabilizer 





-15,16,17, 

right wing 





-18,19,20, 

left King 





-21,22,-23,24 

s^-axis/vertical 






reference 



M22 1 




j 0000 

1 

A+22+1 

4 

-6,2, 1,5 

side A 

0010 

o 

+6 

4 

-8,4, 3,7 

side C 

0100 

4 

+ 11 

4 

-7, 3,2, 6 

side D 

1000 

S 

+ 16 

4 

-5, 1,4, 8 

side B 

1001 

9 

+21 

7 

-6, 2, 1,4,8, -5,1 

side A,B 

1010 

10 

+29 

7 

-5, 1,4, 3, 7, -8,4 

side B,C 

0110 

6 

+37 

7 

“8, 4, 3, 2, 6, -7,0 

side C,D 

0101 

5 

+45 

7 

-7, 3, 2, 1,5, -6, 2 

side D,A 


- 62 . 


TABLE IX. ADDRESS TABLE FOR SQUARE FUSELAGE CROSS-SECTION. 


Code 

1 Address Location Move/Draw List in Core 

0-7 


B ^+1 

B^6 

* *) **) *** **★★) \ 
0 ^ 

B+11 

. B+45 

. . B+37 

0 

8-15 

B+16 

B+21 

B+29 

0 ^ 

0 

0 

0 . 

. .0 


*) A is the address in core of the first location of list 0. 

**) B = A + C where C is the move/draw count of list 0. 

An address value 0 indicates that no move/draw list exists for this 
value of the code. 

**★★) code 11-15 no move/draw lists exist and this part of the table 

might be ommitcd, provided 0 ^ CODE < 10. 






















Value in Units of Half 
the Fuselage Width: d 


Descrip- 

tion 



"p 

"p 


-*8.0 

2.5 

4.0 

hori- 

zontal 

•8.0 

•2 . 5 

4.0 

stabi- 

•5.5 

0. 

4.0 

lizer 

4.0 

HI 

-0.58 


-3.0 


-0.58 


1 

• o 

1 

1.0 

-0.58 

hHI 

•3.0 

1.0 

-0.58 


•3.0 

-6.0 

-0.58 

right 

wing 

4.0 

-1.0 

-0.58 


U. 

0. 


vertical 

0. 

0. 

8.0 

tick- 

marks 

0. 

0. 

•6 .0 

0, 

()• 

•8.0 





































*■' -v 


r* t Cr 
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TABLE XI. MOVE/DRAK CODING LIST FOR PERSPECTIVE VEHICLE SWBOL WITH HEXAGONAL 
FUSELAGE CROSS-SECTION. 


Code 

(Binary) 

000 000 


Code 

(Dec- 

imal) 


i 000 001 

I 000 010 

i 000 100 

; 001 000 
{ 

j 010 000 j 

100 000 1 
— 

! 010 001 

1 

i oil 000 

j 001 010 
! 100 010 
I 100 100 
i 000 101 I 

1 on 001 ! 


i oil 010 

j 

i 

j 101 010 
! 100 110 


100 101 i 


010 101 


List 
Address 
in Core 


M/D 

Count 


A+24+1 


*►125 I 


Coding List 

i 

-9,10,11,12,7,8,9, 

13,14.15, 

-16,17,18,16, 

-19.20,21, 

-22,23,24, 

-25.26,-27,28 

-9.3,2,S 

-12.6,5,11 

-10,4,3,9 

-7.1,6.12 

-S,2,l,7 

-11,5,4,10 

-9, 3, 2, 1,7, -S. 2 
-8, 2, 1,5, 12, -7,1 

-7,1,6,5,11,-12,6 
-12,6,5,4,10,-11,5 
-11,5,4,3,9,-10,4. 
-10,4,3,2,3, -9, 3 


Description 

rear fuselage surface 

vertical stabilizer 

hori zontal stabilizer 

right wing 

left wing 

2 , -axis/vertical 
D 


reference 


side A 
side D 
side F 
side C 
side B 
side n 

side A,B 
side B,C 
side C,D 
side D,E 
side E,F 
side F>A 


-9 , 5 , 2 , 1 » 
-1.7 


7 1 5 ide ,B ,C 


-8, 2, 1,C", 5, 11, -7,1, 1 

- 6,12 

-7,1,6,5,4,10,-12, 

5, -5, 11 

-12,6,5,4,3,9,-11, 
5, -4, 10 

-11,5,4,3,2,8,-10, 

4, -3, 9 

j 

i -10. 4, 3, 2,1,* 


side B,C,D 
side C,D,E 
side D,F-,F 
(side F,F,A 
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of f'O oo^o- t=M. 

.) for codo 43-63 no r.ovo/dtaw 

60 o.,U.«d. provided 0 r C06E 4 43. 











ANALUU 


EAl 560 
linear 


VERTICAL 

DYNAMICS 




EAI 580 
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